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Abstract 

We perform a mean-field study of the binary Bose-Einstein condensate mixtures as a func- 
tion of the mutual repulsive interaction strength. In the phase segregated regime, we find 
that there are two distinct phases: the weakly segregated phase characterized by a 'pene- 
tration depth' and the strongly segregated phase characterized by a healing length. In the 
weakly segregated phase the symmetry of the shape of each condensate will not take that of 
the trap because of the finite surface tension, but its total density profile still does. In the 
strongly segregated phase even the total density profile takes a different symmetry from that 
of the trap because of the mutual exclusion of the condensates. The lower critical condensate- 
atom number to observe the complete phase segregation is discussed. A comparison to recent 
experimental data suggests that the weakly segregated phase has been observed. 
PACS#: 03.75.Fi 

1 Introduction 

Shortly after its first theoretical study [||], the binary mixture of Bose-Einstein condensates (BEC's) 
of alkali-metal atoms in a trap has been realized experimentally 0, and the development of this 
field is now blooming || |J |, ||, [Tj . The binary mixture idea has been extended to trapped boson- 
fermion and fermion-fermion systems Q. This opens the door to studying the rich physics in new 
parameter regimes. There have been several theoretical studies || ||, [7| ; most of them are numerical 
in nature or are from the atomic physics point of view. In the present paper, we perform a mean- 
field-type study based on the nonlinear Schrodinger equation, and obtain a qualitative and in many 
cases a quantitative analytical understanding of a variety of properties of the binary BEC mixtures. 
Along with results from the competition between the healing length and the 'penetration depth', 
as well as from the finite trap effect, we have summarized previously known results in the manner 
of simplified mean-field solutions. In this way, we provide a convenient framework for classifying 
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various excitations in the system, and pave the way for further study of properties of binary BEC's, 
such as the time evolution of the two condensates during a phase segregating process. 
We start from the Hamiltonian formulation of the binary BEC's at zero temperature: 
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Here ipi, with i = 1,2, is the effective wave function of the ith condensate, with the mass rrii 
and the trapping potential C/j. The interaction between the ith condensate atoms is specified by 
Gu, and that between 1 and 2 by Gu- In the present paper all G"s will be taken to be positive. 
The corresponding time independent equation of motion is the well-known nonlinear Schrodinger 
equation ||, obtained here by minimization of the energy, Eq. (1), with fixed condensate atom 
numbers: 
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The Lagrangian multipliers, the chemical potentials fii and /j,2, are determined by the relations 
/ d 3 x\ipi(x)\ 2 = Ni ,i = 1,2, with N{ the number of the ith condensate atoms. 

Experimentally, the trapping potentials {Ui} are simple harmonic in nature. For the sake of 
simplicity and to illustrate the physics we shall consider a square well tapping potential Ui = U: 
zero inside and large (infinite) outside, unless otherwise explicitly specified. 



2 Criteria and Symmetries in Segregated Phases 
2.1 Simplified mean-field solutions 

With the square well trapping potential specified in Sec. 1, the coupled nonlinear Schrodinger 
equations have an obvious homogeneous solution: inside the trap the condensate densities pi = 
|"0i| 2 ; Pi — ' w ^ n V the volume of the square well potential trap, and the chemical potentials 
Hi = Gupi + G12P2 and /i 2 = G22P2 + Gi 2 pi- The corresponding total energy of the system is 
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For a small enough Gu, any variation on top of this solution will increase the system energy. This 
implies that the excitations are stable. Therefore in this parameter regime the homogeneous state 
is the ground state. If the mutual repulsive interaction G12 is strong enough, however, this is no 
longer true. We show below in a mean-field manner that there is an inhomogeneous solution with 
a lower total system energy. 

Let us consider the case of the inhomogeneous state in which the two condensates mutually 
exclude each other. For the moment we ignore the thickness of the interface and the corresponding 
extra energy. In this way we temporarily ignore the derivative terms in Eqs. (2) and (3) in 
determining the effective condensate wave functions. We call this situation the simplified mean- 
field approach, which is a useful one that has already given us a lot of physical insights QXJ . Let V* 
be the volume inside the trap occupied by the condensate i. We have \ipi\ 2 = p iQ = Ni/Vi and the 

1 N 2 

total energy of the inhomogeneous state E in = ~J2 i=12 Gii^-. Minimizing E in with respect to V\ 
or V 2 under the constraint V% + V 2 — V, we obtain the spatial volume occupied by each condensate: 



1 1 / G22 N2 



-V , Vo 



1 



-V . 



The corresponding condensate densities are 
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and the chemical potentials /ij = Ga p i0 . We note here /iipio = P2P20- The total energy for this 
inhomogeneous state is 
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The energy difference from the homogeneous state is then 

AE = E in — E ho = — Gu — v G11G22 
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This equation reveals that for a large enough mutual repulsive interaction, that is, if 
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the inhomogeneous state has a lower total energy. Hence, the inhomogeneous state, the phase segre- 
gation state, will be favored for a large mutual repulsive interaction Gi2- We note that this criterion 
for the mutual repulsive interaction strength is independent of the condensate-atom numbers as well 
as of the trap size. We shall return to this point below. The critical value for G 12 , Eq. (7), has 
been found using a stability analysis from the excitation spectrum in the homogeneous state ||, 
while it is obtained here from a simple energetic consideration. 
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2.2 Interface profile 

In the absence of the derivative terms in the determination of the condensate wave functions studied 
above, the shape of the boundary between the two condensates in the phase segregated state will 
take any form. Now, we consider the effect of the derivative terms. Their inclusion will make the 
thickness of the interface finite and introduce a finite interface energy, the surface tension. The 
presence of the surface tension will fix the shape of the interface between the two condensates by 
minimization of the total surface energy. First, we look for the condensate profiles at the interface 
with a finite thickness. We rescale the effective wave functions by the values deep inside their own 
condensates: 

A = fiy/PiO ■ 

We shall assume here that condensate 1 occupies the region in the trap specified by z > 0, and 
condensate 2 occupies the region specified by z < 0. The interface plane is z = 0. 

Deep inside the region of condensate 1, we have f\ — > 1 and fa — > 0. In this region, using 
Pi = GuPi and only keeping the leading contributions, the coupled nonlinear Schrodinger equations, 
Eqs. (2) and (3), becomes 
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Here fa — 1 + 6 fa. These equations may be written in the following more suggestive form 
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From Eq. (9) we find the density profile of condensate 2 in condensate 1 is 

fa(z) = / 2 (0)<r£ , (11) 

with the 'penetration depth' 

A -' = i m • (12) 

which is the length scale for condensate 2 penetrating into 1. Similarly, for the penetration of 
condensate 1 into condensate 2 in region z < 0, we have fa(z) = fa(0)e z ^ Al , with the penetration 
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depth Ai = ii/ \J jg^g 22 ~~ 1- The healing length scale £ here describes the ability of a condensate 
to recover from a disturbance, similar to the same length scale in superfluid helium 4. The newly 
introduced length scale here, the penetration depth A, describes the degree of the mixing between 
the two condensates. Obviously, as G12 — > y/GuG^, the penetration depth goes to infinity, in 
coincidence with the disappearing of the phase segregation. 

It is also useful to study the behavior of condensate 1 in the region z > 0. Deep inside condensate 

1, 5fi is small. If /| approaches zero faster than 5fi, that is, /f < 5fi, the last term in Eq. (8) 
may be dropped, and we have 

5h = 5/i(0) e-^ 1 . (13) 

Here it is the healing length £1 of condensate 1, not the penetration depth A2, that determines the 
profile of condensate 1. The validity of self-consistency for this solution requires A2 < \/2£i. In this 
parameter regime the mutual repulsive interaction is so strong that condensate 1 stays away from 

2, which is similar to the Meissner state where the magnetic field is completely excluded outside of 
a bulk superconductor. We shall call this parameter regime the strongly segregated phase. In the 
opposite limit, that is, A 2 > \/2£i, 5fi will be determined by f2 through Eq. (8): 

2 _ iS^J VG11G22 

The only relevant length scale here is the penetration depth, which is larger than the healing length. 
There is still a considerable mixing of the two condensates in this parameter regime, which we shall 
call the weakly segregated phase. 



2.3 Surface tension 

With the inclusion of the gradient terms in the determination of the condensate wave functions, the 
presence of the interface will cost a finite amount of energy. The surface energy per unit area, the 
surface tension, may be defined as a = AE S /S. Here S is the interface area. The energy difference 
AE S may be calculated in the following manner: In the presence of the interface at z = we first 
solve the full Eqs. (2) and (3) with derivative terms for the condensate wave functions fy, calculate 
the corresponding total energy from Eq. (1), then subtract from this total energy by the amount 
given by Eq. (5), the total energy of the system in the segregated phase without the effect of the 
derivative terms in Eqs. (2) and (3). Specifically, the energy difference AE S is 
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Because of the normalization condition / d 3 x Pioff = iVj, from Eqs. (2) and (3), we obtain the 
surface tension as 
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Though Eq. (15) is a general expression for the surface tension, to gain a concrete understanding, 
we consider the case that the two condensates have an identical set of parameters: Aj = A and 
£i = £. In the strongly segregated phase of £ >> A/y/2, Eq. (15) gives 



o- = -j= \JG U G 22 P10P20 , (16) 

which is independent of the penetration depth and the mutual repulsive interaction. In the weakly 
segreg ated phase of £ « A/V2, Eq. (15) gives 




which goes to zero as A — > 00. We note that this occurs when G\% — > y/GuG^, in agreement with 
our above mean-field analysis of the phase segregation. The existence of the finite surface tension 
leads to another branch of gapless excitations, the interface or surface mode, which we are not going 
to discuss here. 



2.4 Finite trap size effect I: broken symmetry ground state 

Now we consider the effects of finite surface tension and the finite trap size. For a very large system, 
it is known that the minimization of surface energy leads to the minimum surface area, whose shape 
is usually spherical in three dimensions (3D) and circular in 2D. For a finite size trap, one might 
expect that the shape of the ground state of the binary BEC mixture should take the same symmetry 
of the trap, particularly if it is cylindrically or spherically symmetric. We show here that this may 
not be true in a segregated phase, and the condensates may break the cylindrical symmetry defined 
by the trap. To demonstrate the essential physics, we consider the case of two condensates in two 
spatial dimensions with an identical set of parameters: Aj = A, £3 = £, and pio = po- Let R be the 
radius of the square well potential trap. Supposing that condensate 1 occupies a circular area with 
a radius R! = ^j=; the associated surface energy is the length of the interface 2nR' times the surface 
tension a : 

E s = V2 nRa . 

On the other hand, supposing that each condensate occupies a half circular shape of the trap, the 
corresponding surface energy is 

E b = 2Ra , 

which is much lower than the circular shape with the same symmetry as the trap, because the 
interface length here is shorter than that with the circular symmetry. Therefore, each condensate 
will take a different shape than that of the symmetry of the trap. The broken cylindrical symmetry 
state occurs, discovered first numerically 0. 

Though the circular symmetry is broken for each condensate, in the weakly segregated phase 
where A > v^2£, the total density p\ + P2 still appears circularly symmetric, and it retains the 
symmetry of the trap. In the strongly segregated phase where A < \/2£, the two condensates tend 
to avoid each other. In this regime the circular symmetry of the total density profile is broken. 
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Here we give a heuristic demonstration of such different behaviors of the total density profile in two 
segregated phases. We first show that there is a ditch in the total density profile at the interface in 
the strongly segregated phase. Following the analysis above, we take z = as the interface position. 
Condensate 1 (2) occupies the z > (z < 0) region of the trap. In the region z > 0, the condensate 
densities take the forms pi(z) = p (l — b\e~^ z ^) and p 2 {z) = po b 2 e~ 2z / A , consistent with Eqs. (8) 
and (9). Similarly, in the region z < 0, Pi(z) = p b 2 e 2z ^ A and p 2 (z) = p (l — b 2 e^ 2 ~ z ^). Here b\ and 
b 2 are two numerical constants. To determine b\ and b 2 , we make use of the fact that because Eqs. 
(2) and (3) are second order differential equations, the solutions and their first order derivatives 
must be continuous. This immediately gives us two algebraic equations at the interface z = 0: 

1 - 6i = b 2 , v^fei/e = 26 2 /A . 

The numerical constants are then 6i = y/2Z/(y/2£ + A) and b 2 = A/(v / 2^ + A). Evidently, the total 
density pi(z) + p 2 {z) has its minimum value po 2A/(y / 2£ + A) at the interface z — 0, which goes 
to zero as A — > or, equivalently, as G 12 — > oo. The total density has a ditch at the interface, an 
indication of broken symmetry in the strongly segregated phase. For the weakly segregated phase, 
replacing \/2^ by 2 A and following the same procedure we obtain bi = b 2 = 1/2 and find that 
the total density remains constant in the trap, not affected by the phase segregation. There is no 
broken symmetry for the total density in the weakly segregated phase. 

We note that in terms of interactions the weakly segregated phase is specified by 1 < G 12 /y/G n G 22 < 
(1 + l/y/2), and the strongly segregated phase by Gi 2 /\/G n G 22 > (1 + \/y/2). 

Another interesting feature of the finite surface tension is the floating of the condensate droplet. 
This may be regarded as a special case of the symmetry broken state. Suppose one condensate, 
say condensate 2, has a particle number much smaller than that of condensate 1, but still large 
enough to have a well-defined interface and surface. Then condensate 2 may form a droplet inside 
condensate 1 in a phase segregated regime. Because of the finite surface tension, this droplet may 
move to and stay at the edge of the trap to reduce the common boundary length in 2D (or area in 
3D) between condensates 1 and 2, to minimize the surface energy. This tendency may be called the 
floating of the condensate 2 droplet. 

2.5 Finite trap size effect II: lower critical condensate-atom number 
for phase segregation 

In the above analysis, we have implicitly assumed that the thickness of the interface is much smaller 
than the trap size, such that we have a well-developed phase segregation. We examine the limitation 
of this assumption here. According to the above analysis, the penetration of condensate 2 into 1 is 
determined by A, and the recovery of condensate 1 from the presence of the interface by £ or A in 
the strongly or weakly segregated phase. As an estimation, we may have the interface thickness l 12 
as 

/i2~£ + A. (18) 

This implies that l\ 2 oc 1/ v /pIo. When the interface thickness is larger than the trap size, l\ 2 > R, 
we do not have a well defined phase segregated state. A conclusion inferred from this is that, if the 
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condensate atom number is too small, there is no complete phase segregation in the trap even with 
a strong mutual repulsive interaction between the condensates. 

We should point out that the analyzes in Sees. 2.4 and 2.5 are based on the assumption that 
the two condensates have the same sets of healing length and penetration depth. This implies that 
the surface energy contributions caused by the edge of the square well trap are the same for both 
condensates. Hence we have ignored their effect in the determination of the symmetry of the ground 
state. In general, this may not be true, and one has to consider the competition among four lengths: 
two healing lengths and two penetration depths, as well as the surface energy from the trap edge. 
This will generate an even richer physics than what has been presented above. We believe our 
above discussions have provided the framework for further explorations. One such example will be 
discussed in Sec. 3. 



3 Numerics and discussions 

In terms of the atomic scattering lengths of condensate atoms an, the interactions are Ga = 
ATih 2 aa/mi. The typical value of an for 87 Rb is about 50 A. The typical density realized for the bi- 
nary BEC mixture is about pio ~ 10 14 /cm 3 . Hence the healing length is £ = J (ft 2 /2m)(l/GuPio) — 

^Jl/(87iauPio) ~ 3000 A. For the different hyperfme states of 87 Rb, we may take Gi 2 / VG11G22 = 
1.04, whose precise value is uncertain and may be smaller. Then the penetration depth is A = 

i/ \J G12/ \fGi\G22 — 1 ~ 1.5/iHi. The length scale for the harmonic trap potential determined by 
oscillator frequencies ranges from 1.3 to 4 /im, and the condensate occupies a region with a diam- 
eter of about 20 pm for about 10 6 atoms, which is comparable to or larger than the size of the 
interface thickness. Hence it is reasonable in practice to apply the mean-field results obtained for a 
square well potential to the case of a harmonic trapping potential, because on the scale of interface 
thickness the trapping potential appears smooth. The measurement of the penetration depth can 
be used to determine the mutual repulsive interaction in the 87 Rb system, which is in the weakly 
segregated phase, because Gu is believed @ to be slightly larger than y/G\\G 22- 

Now we consider whether or not the ground state symmetry of the hyperfme states 87 Rb can be 
broken. As pointed out at the end of Sec. 2.5, there is a complication arising from the trap edge 
surface energy contribution, due to the difference in healing lengths or interactions. In accordance 
with the experimental situation and for the sake of simplicity of analysis, we take two condensates 
having an equal number of particles. The differences in interactions are small: {Si = (Ga — 
G\2)/Gi2, i = 1,2} are close to zero, and their sum Si + 62 ~ 0. A condensate near the edge of the 
square well potential is identical to the strong segregated phase, because there is no penetration 
into the hard wall. The surface tension at the trap edge can be readily evaluated according to an 
expression similar to Eq. (15) in the strongly segregated phase but with only one condensate. The 
difference in those surface tensions is then, up to the first order in Si, 

Act = cti - cr 2 = — x \(S 2 - Si) . 



Here we have used Eq. (5) that V G22 P20 — \fG~11 Pio- If there were no symmetry breaking in the 
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ground state, the condensate with the lower surface energy at the trap edge, say condensate 2, will 
stay close to the trap edge. Condensate 1 stays inside. If there is a symmetry breaking such that 
each occupies half of the trap as discussed in Sec. 2.4, then condensate 1 will get in contact with 
half of the trap edge. There is an increase in total energy as 

AE edge = TrRAa = ^R^Gnplo X (fc - <*i) • (19) 

For the symmetry breaking to occur, this surface energy cost from the trap edge must be smaller 
than the interface energy gained, which is 

AE inter f ace = E s — E b = ^—R£iGiipl x (82 — 81) . (20) 

Here we have used £j/A$ = (v2/4) (5 2 — 5i) and 8\ + 5 2 = 0. One can readily check now that indeed 
AE inter face > AE e d ge - We note that in the case of the harmonic trapping potential, the density of 
a condensate is smaller near the edge, which gives a even smaller edge surface energy contribution. 
Therefore the symmetry breaking will occur in the system of hyperfine states of 87 Rb according to 
the present analysis. We believe this is precisely what was observed in a very recent experiment [BJ. 
Furthermore, since it is in the weakly segregated phase, one should also expect that the symmetry 
break occurs only for the density of each individual condensate, not for the total density. Again, 
this is what was reported in Ref. 
onlinecitefurther. 

Given the size of the trap to be on the order of 20 jum, the interface thickness must be smaller 
than this length to have a well-defined segregated phase. For the different hyperfine states of 
87 Rb, this implies that, according to Eq. (18), a lower critical number N c of condensate atoms 
N c = l/(Gi2v / G'iiG'22 — 1) V 1 / 3 /(%-ncLu) ~ 4000. We should point out that the precise value of the 
lower critical atom number depends on details of a realistic trapping potential, such as the oscillator 
frequency and the anisotropy ratio. For a condensate-atom number of less than this value, there is 
no complete phase segregation. 



4 Conclusion 

From the mean-field analysis, by tuning the strength of the mutual repulsive interaction we have 
found that there are two segregated phases: The interface profile is determined by the penetration 
depth in the weakly segregated phase , and by the healing length in the strongly segregated phase. 
The broken cylindrical symmetry state starts to appear in the weakly segregated phase for each 
condensate, and persists into the strongly segregated phase. For the total condensate density, the 
cylindrical symmetry is maintained in the weakly segregated phase, and disappears in the strongly 
segregated phase. We have also found that a condensate droplet inside another condensate in a 
segregated phase tends to move to and to stay at the trap boundary, and that if the condensate atom 
number is smaller than a critical value, there is no well-developed phase segregation. A comparison 
between the present results and a recent experiment has suggested the weakly segregated phase has 
been observed in 87 Rb condensates. 
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